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SUMMARY 

A published model and ancillary computer programs were used for a 
comparison of theory with an actual conical, pressurized-fluidized-bed 
combustor burning caking bituminous coal and using limestone to reduce 
sulfur dioxide emission. It was necessary to make a number of changes 
to the source programs to get meaningful results. In addition, the 
inputs of limestone density, adjustable elutriation parameter; and heat 
transfer coefficient were adjusted to compensate for limitations in the 
model. Theoretical bed pressure drop was in good agreement with 
experiment. The burnable carbon elutriated indicated that, contrary to 
the assumption in the theory, the exhaust port was below the transport 
disengaging height, resulting in disagreement between theory and 
experiment. The observed nitrogen oxides emission rate was about half 
the theoretical value. There was order-of-magnitude agreement on 
sulfur dioxide emission rates. Recommendations for improving the model 
are given. 

INTRODUCTION 

The pressurized-fluidized-bed coal combustor is being investigated 
by the Department of Energy, the utility industry, and several 
laboratories with the ultimate purpose of achieving clean coal 
combustion in high-efficiency central -station powerplants (refs. 1, 

2, and 3). In such systems, the ability to scale up small combustors 
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to larger sizes is of great importance. Ideally, the best way to scale 
up such combustors would be with a realistic physical model as opposed 
to multiple regression correlations or empirical scaling laws. The 
primary purpose of this report is to give preliminary comparisons 
between a physical model (ref. 4) and a conical pressurized-fluidized- 
bed combustor (refs. 5 and 6) and attempt to explain any 
discrepancies. A secondary purpose of this report is to provide 
additional information for anyone using the two complementary computer 
programs published in reference 4. The scope of this work did not 
include a systematic check of the equations used in the model nor the 
computer programs. A number of discrepancies are reported which were 
discovered in the course of attempting to run the programs. Most of 
these were corrected. 

Previous work modeling fluidized bed combustors has been somewhat 
limited. Horio and Wen (ref. 7) used a modified bubble assemblage 
model with distributions of coal particle size and limestone conversion 
to calculate sulfur dioxide removal efficiency by limestone additive. 
Horio and Wen (ref. 8) later published a treatment of elutriation of 
fines, solids mixing, combustion efficiency, and bed temperature 
profile. Chen and Saxena (ref. 9) included both sulfur dioxide removal 
efficiency and combustion efficiency and were first with a three-phase 
comprehensive model. Concurrent work in modeling includes that of the 
MIT Energy Laboratory staff (ref. 10) and that of Rajan, Krishnan, and 
Wen (ref. 11). 
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All the previous modeling work mentioned applies to fluidized bed 
combustors of constant cross-sectional area. This report is the first 
comparison of a variable cross-sectional area (conical) fluidized bed 
combustor with a model (ref. 4). The conical bed has the advantage 
that it minimizes elutriation (ref. 5). The realism of the model was 
improved by using the actual bed particle density, which was measured 
three ways. 

APPARATUS 

Conical, Pressurized, Fluidized Bed Combustor 
The pressurized fluidized bed combustor and associated systems are 
shown in a simplified schematic in figure 1. Geometrically, the 
combustor consisted of a lower 24.7 inch high cylindrical section of 
8.8 inch inside diameter, a middle 7.3 inch high truncated-cone section 
of 10.48° cone half angle, and an upper 79.9 inch high truncated cone 
section of 3.40° cone half angle, bringing the inside diameter at the 
top to 21 inches. There were six ports spaced vertically on the side 
of the combustor. A solids removal auger could be located in any one 
of the ports to maintain the level of the bed inside the combustor no 
higher than that port. 

A mixture of coal and limestone was injected into the bottom of the 
combustor using high pressure air as a transport media. The fuel (coal 
plus limestone) flow was controlled by the rotational speed of the fuel 
metering screw, which was calibrated in terms of flow. Pressurized air 
at ambient temperature flowed into the bottom of the reactor through a 
distributor containing a total of 36 holes 1/8 inch in diameter in nine 
bubble caps. 
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Heat was removed from the bed via banks of water-cooled horizontal 
heat-exchanger tubes. Three banks of seven tubes per bank were 
installed in the bottom two feet of the bed. 

The limestone-to-coal ratio was determined from the rotation rates 
of the coal and limestone metering screws, which were calibrated for 
flow rates. The air feed rate was determined with venturi meters. 

Other details are given in reference 5. 

Cylindrical Cold-Flow Atmospheric Fluidized Bed 

A cylindrical fluidized bed was used to measure the minimum 
fluidization velocity and void fraction of bed material taken from the 
conical combustor. This was part of one method of determining particle 
density of spent bed material. The apparatus is shown in figure 2. 

The fluidizing gas was dry air. The inlet temperature of the air was 
measured with an alcohol-in-glass thermometer in an oil-filled 
thermometer well. Flow was manually controlled with a throttle valve 
and read on a rotometer. The reading was corrected for pressure 
variations by means of a bourdon-tube pressure gage upstream. The bed 
was contained in a transparent plastic column with a fine wire mesh 
distributor near the bottom. The top of the bed was at ambient 
pressure. Bed pressure drop was measured with a diaphragm pressure 
gage. 

MODEL 

The theoretical model for the conical, pressurized fluidized bed 
coal combustor is given in reference 4 (the "Level II" model), but will 
be reviewed briefly here. 



5 


Hydrodynamics 

It is assumed that "fast" bubbles are present (rising velocity of 
bubble greater than gas velocity in emulsion). Hence, the bubbles have 
clouds. The bubble size is given by the correlation of Mori and Wen 
(ref. 12) modified for a combustor of varying cross-sectional area. 

Char mixing is calculated from a single-phase back-flow multicell 
model. Limestone is assumed to be completely mixed. Gas flow is based 
on a two-phase theory, where the emulsion is one phase and the bubbles 
and clouds are the other phase. A gas interchange coefficient between 
the two phases is included. 

Coal Combustion 

A spherical particle model is employed for coal combustion. Only 
the lean case is treated. The hydrogen and oxygen volatile immediately 
upon injection of coal into the combustor, not changing the diameter of 
the resultant char. The diameter is gradually reduced by burning with 
oxygen, with the rate determined by the surface rate of chemical 
reaction and gas diffusion. Ash particles break off as the char 
burns. Carbon, nitrogen, and sulfur in the char are assumed to be 
released or used at the same rate as the char burns. 

Limestone Chemical Kinetics 

A porous constant diameter sphere model is employed for limestone 
chemical kinetics. The limestone is assumed to calcine as soon as it 
is injected into the combustor, producing pores. The reaction rate 
with sulfur dioxide is given by Borgwardt's (ref. 13) model. 
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Nitrogen Oxides Kinetics 

Nitrogen oxides are assumed to be produced only from coal-bound 
nitrogen. A numerical fit to the experimental results of Ruth (ref. 

14) was made, assuming the quantity of nitrogen oxides produced is 
directly proportional to the BTU's of the coal actually burned and 
nonlinearly proportional to the excess air. 

Heat Balance 

A multi cell single-phase backmix model is employed. The overall 
heat transfer coefficient to the heat exchanger tubes must be given. 

The heat transferred to the walls and radiated from the top of the bed 
is assumed to be zero. 

Elutriation 

Elutriation of char, ash, and limestone are treated differently. 

For char it is assumed that the combustor exhaust gas port is above the 
transport disengaging height, and one of three empirical correlations 
may be selected for the elutriation rate. The fraction of ash 
elutriated is not calculated so it must be given as an input. The 
limestone is assumed not to elutriate. 

COMPUTER PROGRAMS 
General Description 

The detailed model (Level II) is too complex even for a high-speed 
digital computer. Hence a simplified model is used to estimate the 
loss of burnable carbon due to elutriation and the combustion 
efficiency. The program for this simplified model is called Level I. 
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Level I Program 

This program assumes plug flow or back mixing of gas (plug flow was 
assumed in this report). It assumes complete mixing of solids and an 
isothermal cylindrical bed. The coal fed is assumed to have a size 
distribution. The bed temperature must be given. The size 
distribution of the char is solved for, taking into account combustion, 
elutriation, and solids withdrawal. This requires three nested 
iteration loops. The inner loop is for (symbols are given in 
Appendix A), which is related to the ratio of coal particles fed to 
char particles withdrawn. The middle loop is for(p, which is the mean 
residence time of limestone. The outer loop is for?/^, which is the 
combustion efficiency. Outputs include the loss of burnable carbon due 
to elutriation. 

Level II Program 

This program assumes a bed of any geometry, size distribution of 
coal, and size distribution of limestone. Inputs include 7^^ and the 
loss of burnable carbon due to elutriation, both from Level I. It 
solves the rest of the equations described in the MODEL section of this 
report, including bed temperature and burnable carbon concentration as 
functions of height. 

Changes to Programs 

A few changes to the two programs were made merely for convenience, 
such as changes to output format. However, other changes were 
programmed that affected the outputs. These are given below, where all 
card numbers refer to the numbers on the left on pages 62-91 of 
reference 4. 
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Level I Program 

Changes are given in Table I. Here MAIN refers to the main Level I 
program. These changes were not 100 percent successful. In 45 cases, 

did not converge 348 times, but always converged on the final 
iteration. This is believed related to the change in subroutine POP to 
eliminate incorrect negative values of (j>. The differential equation 
for is solved from y = 0 to y = 1. If it were solved from y = 1 to 
y = 0, it is believed negative values of would never occur even 
without the change in POP tn Table I. Also, in two cases the 
exponential in POP was out of range. Attempts to correct this were 
abandoned. 

Level II Program 

Changes are given in Table II. These changes appear to have been 
successful . 

PROGRAM TEST CASES AND ADJUSTMENT OF CONSTANTS 
Test Cases 

Level I Program 

A test case is given in reference 4, but the input contains some 
errors: SCF should read XCF, DPF(7) should be 0.1003 instead of 0.1, 

and EXAIR should be 0.06 instead of 0.13. With these changes (but not 
the changes in Table I), the outputs from an IBM-360 computer at NASA 
checked reference 4 (see Table III) although the residence time 6^ did 
not converge. With the changes in Table I, the outputs were different 
(see Table III), although the only outputs used as inputs in Level II 
(ETC (combustion efficiency) and ELOSS (carbon elutriated)) were not 
much different. 
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Level II Program 

A test case is given in reference 4, but the input contains some 
errors: ATB(l) and ATB(2) should be 4180.6 instead of 4181.0, ZB(4) 

should be 150. instead of 200.0, ZHE(4) should be 300. instead of 
200.0, ZDIS(l) should be 66. instead of 0.0, DNZL should be .37338 
instead of 0.3734, DEAV should read DZAV, and EXAIR should be .174 
instead of 17.4. With these changes (but not the changes in Table II) 
the outputs in reference 4 were obtained. With the changes in Table 
II, the output was as given in Table IV. The only large change is in 
ANOX(NO^). 

Preliminary Results and Comparisons 

For comparison, five experimental steady-state tests with a total 
duration of 20 hours and 7 minutes were run under the same conditions 
with the discharge solids removal auger in port #3 ^see figure 1). 
Caking bituminous coal from the Pittsburgh #8 seam was used with Grove 
City, Virginia limestone. The material loaded in the bed prior to the 
tests was bed material from previous tests using the same kind of coal 
and limestone. The results of the tests were averaged and are given in 
Tables V and VI as well as figure 3, figure 4, and Appendix B. The 
methods of calculating or obtaining inputs for Level I and Level II 
programs are given in Appendix B for the quantities which were not 
changed in later calculations for the five averaged tests. 
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The input quantities which were changed later were HLMF (bed 
height), DPI (limestone diameter), BETA (adjustable elutriation 
parameter), DB (bubble diameter), and RHOLS (limestone density) in 
Level I; and ELLOSS (carbon elutriated), RHOAD (limestone density), 

ETCA (estimated combustion efficiency), NDPAD (number of limestone 
sizes), DPADF (limestone size), FRACTA (fraction of limestone of a 
size), and UHEAV (heat transfer coefficient) in Level II. We shall 
call the preliminary cases for the five averaged tests Case la. Case 
lb, and Case Ic for elutriation correlations from references 15, 16, 
and 17, respectively (corresponding to lELUTR of 1, 2, and 3 in Level 
I). For Cases la, lb, and Ic, the above input quantities were 
calculated or obtained as given in Appendix C. 

Outputs from Level I and II programs are given in Tables V and VI, 
respectively, along with observed values (observed ELOSS (carbon 
elutriated) was obtained using a fly ash burnable carbon fraction 
determined by measuring carbon on a sample previously washed with 
dilute phosphoric acid to eliminate all carbonates). The calculated 
fraction of burnable carbon elutriated (ELOSS) was far lower than 
observed for all three elutriation correlations (refs. 15-17). 

However, Case la (lELUTR = 1 so reference 15 used for elutriation) was 
closest (see figure 3). 

For Case la, the bed temperature was too high, sulfur retention was 
too low, fractional conversion of additive was too low, and sulfur 
dioxide and nitrogen oxides in the exhaust were too high compared to 
observed values. In addition, there was slight agreement between 
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calculated and observed size distribution of elutriated carbon (see 
figure 4), the observed carbon being much smaller than the predicted 
carbon. 

Changes to Input Data and Comparisons 

Case 2a 

The particle density of a bed of porous particles such as limestone 
in various degrees of calcination and sulfation is difficult to measure 
because it is supposed to be based on the actual mass of the particles 
and a volume that includes both the solids and the pores but not the 
interstices. Different fluids used in measuring the particle density 
go into the pores to various degrees. To find a realistic particle 
density for the spent bed of the conical fluidized bed combustor, the 
measurement on the spent bed was made in three ways: (1) a pyknometer 

with Jet A kerosene gave a density of 2.73 g/cm ; (2) a mercury 

3 

porosimeter gave densities of 2.29 and 3.05 g/cm at pressures of 1.8 
and 15000 psia, respectively; (3) a cylindrical cold flow fluidized bed 
(figure 2) gave a minimum fluidization velocity which in conjunction 

3 

with a dry sieving analysis gave a particle density of 2.62 g/cm . 

The last of the three results (2.62) was then compared with the 
particle densities calculated for the bed (calcined limestone) in Level 
I and II programs. The calculated particle densities were too low. 
Consequently, the fictitious limestone densities needed to give a bed 

3 

particle density of 2.62 g/cm were calculated and were 4.89 and 4.57 
g/cm for Level I and II, respectively (see Table V and VI). In 
Level I, the spend bed dry sieving analysis and equation (Cl) gave a 
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DPL (limestone size) of 0.0793 cm. From equations (C7) - (Cll), DB 
(bubble diameter) was 13.75 cm. From equations (C12) - (C13) HLMF (bed 
depth) was 62.4 cm. In Level II DPADF (limestone size), NDPAD (number 
of limestone sizes), and FRACTA (fraction in limestone size) were taken 
from the dry sieving analysis of the spent bed (columns 1 and 4 of 
Table VII). The resulting Level I and II outputs are given in Tables V 
and VI, respectively. ETS (sulfur capture), FS (limestone conversion), 
and XG0(3) (SO^ in offgas) were much closer to the observed values, 
while ELOSS (carbon elutriated), ETC (combustion efficiency), and T 
(temperature) were slightly worse. 

Case 3a 

To get Level I to reproduce the observed ELOSS (carbon elutriated), 
BETA (elutriation parameter) was adjusted to 25. The Level I and II 
outputs are given in Tables V and VI, respectively. 

Case 4a 

To get T (temperature) at 36 cm above the top of the distributor to 
agree with the observed value, UHEAV (heat transfer coefficient) was 
adjusted to 1.083x10 cal /cm sec K. This compensated for the 
lack of wall heat transfer and radiant heat transfer from the top of 
the bed in Level II. The Level I and II outputs are given in Tables V 
and VI, respectively. 

The constants for Case 4a were used for Cases 5 to 30 below, but in 
Level I the variables WCOAL (coal flow), CABS (calcium to sulfur 
ratio), EXAIR (excess air), P (pressure), TK (temperature), OT (bed 
diameter), HLMF (bed depth), DB (bubble diameter), and PASH (ash 


/ 
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elutriated) were changed to conform to the results of the new 
experiments described below. In Level II, ELLOSS (carbon elutriated), 
ETCA (estimated combustion efficiency), HLF (bed height), PAV (bed 
pressure), TAV (temperature), WCOAL (coal flow), CABS (calcium to 
sulfur ratio), PF (pressure below distributor), EXAIR (excess air), and 
GZCO (ratio of CO to C in coal) were also changed to conform. 

END RESULTS AND DISCUSSION 

Twenty-six new tests of about four hours each were run at different 
conditions but with coal and limestone from the same source. The 
calcium-to-sulfur ratio, excess air, bed pressure, bed depth, and bed 
temperature were varied. Level I and II programs were used to predict 
observable quantities for the 26 tests (Cases 5 to 30). The results 
are given in figures 5-12. 

Figure 5 gives observed and calculated bed pressure drop. The 
agreement tends to justify the use of the fictitious limestone 
densities. The wild point is probably due to an instrumentation 
problem. 

Figure 6 gives observed and calculated fraction of burnable carbon 
elutriated. The scatter is great. The four most likely reasons for 
the scatter are: (1) the exhaust port is below the transport 
disengaging height so it gets splashed by particles from bursting 
bubbles. This is confirmed by figure 13; (2) char attrition was 
neglected; (3) Size segregation with height due to hydrodynamic effects 
was neglected (this is very pronounced for spent bed material in the 
cylindrical bed, figure 2); and (4) char thermal decrepitation was 
neglected. 
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Figure 7 gives observed and calculated combustion inefficiency. 

Most of the combustion inefficiency is due to solid carbon elutriated 
or discharged. The observed values are based on bomb calorimeter data 
for both fly ash and discharge. The scatter is partly due to all the 
observed burnable carbon fractions in the discharge being measured as 
zero or less due to basic problems with using a bomb calorimeter for 
this application (negative carbon fractions were arbitrarily set equal 
to zero). 

Figure 8 gives observed and calculated nitrogen oxides emission 
rate. The calculated values tend to be too high. Ruth (ref. 14) 
states that "Larger combustors can be expected to emit even less NO^ 
because of more uniform temperatures and lower combustion 
intensities." Since Ruth had a 4.5 inch i.d. combustor compared to our 
8.8 inch i.d. combustor, this is consistent with figure 8. In 
addition, Beer et al. (ref. 18) have pointed out that nitric oxide is 
partially reduced by reaction with char in the bed. This reduction 
would be more pronounced in our combustor than Ruth's because on the 
average our bed was deeper and, due to the conical shape, had a lower 
superficial velocity near the top. 

Figure 9 shows observed and calculated bed temperature. The reason 
for the scatter is not known, but may be due to difficulties in 
measuring the fuel feed rate by means of the fuel metering screw. 



15 


Figures 10, 11, and 12 give observed and calculated sulfur emission 
rate, efficiency of sulfur capture, and calcium utilization, 
respectively. These are all related. The agreement with theory was 
very limited. Multiple regression analyses did better, showing 
moderate scatter for sulfur emission rate and efficiency of sulfur 
capture about the regression equations, but large scatter for calcium 
utilization. The extremes of sulfur emission rate are so wide that one 
or both of two mechanisms are suspected: (1) raw coal is transported 

by bubbles to the top surface of the bed once in a while, or (2) hot 
spots are formed in the bed, causing the limestone to lose SO^ there 
(ref. 19). The scatter in calcium utilization is probably due to the 
duration of each test being considerably shorter than the residence 
time of the limestone (24 to 174 hours) so the condition of the 
limestone is influenced by previous tests. 

In many cases the bed temperature from the Level II program 
increased with height all the way from the 36 cm level to the top of 
the bed. The increase was as much as 8.8%. Experimentally the 
increase never exceeded 0.5%, and generally the temperature decreased 
with height or was essentially constant. The difficulty is there is no 
wall heat transfer in the Level II model (and radiant heat transfer 
from the top of the bed was neglected but may be important). 



16 


Some of the lack of agreement between Level II output and observed 
values may be due to the impossibility of making a valid adjustment to 
the wake-to-bubble volume ratio f^ in the program due to insufficient 
experimental data. The quantity f determines the degree of solids 

W 

mixing and could be adjusted if burnable carbon as a function of height 
were measured, or wall heat transfer was included in the program. 

Another contributor to lack of agreement between Level I and II 
output and observed values may be the bubble sizes assumed in the 
programs, which, as the authors point out (ref. 4), have never been 
verified for conical beds with heat transfer tubes in the bed. It 
appeared that the calculated bubble sizes were larger than observed 
with closed-circuit television (fig. 1). 

In general, the results from Level I and II were not as good as 
multiple regression analyses in correlating the experimental data, 
although one would hesitate to trust a multiple regression analysis for 
gross extrapolations. 

CONCLUSIONS AND RECOMMENDATIONS 

1. Bed pressure drop from theory and experiment were in good agreement. 

2. Measured burnable carbon elutriated was 21 times the theoretical 
value until the value of an adjustable parameter was changed from 1 to 
25. 

3. The observed average size of char elutriated was about one-third 
the theoretical value. 
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4. The observed nitrogen oxides emission rate was about half the 
theoretical value. 

5. There was order- of-magnitude agreement between theory and 
experiment on sulfur dioxide emission rate. 

6. The published digital computer program is not fully operational, 
and the computer results were inferior to multiple regression analysis 
for correlating the data. 

7. The theoretical burnable carbon elutriated should be improved by 

including a correlation for exhaust port below the transport 
disengaging height and perhaps one or more of the following: (1) char 

attrition, (2) char thermal decrepitation, or (3) char segregation by 
hydrodynamic forces. 

8. The bubble size correlation should be verified for conical beds 
with heat exchanger tubes in the bed. 

9. The theory should include wall heat transfer (and perhaps radiant 
heat transfer from the top of the bed). 

10. The theoretical mixing should be adjusted by means of an empirical 
factor, preferably by measurement of burnable carbon as a function of 
height in the bed. 
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APPENDIX A - SYMBOLS 

Mathematical Level I Level II 

FORTRAN FORTRAN 

At Area of top of bed 

ALAMA Average value of (number 

basis) 

ALAMV Average value of (number 

basis) 




AHE 

Specific heat transfer area of 
heat exchanger tubes 



ANOX 

Mole fraction NO^ in the 
effluent gas 



ATB 

Cross sectional area of 
combustor at specified 
heights 

^cf 

BCl 


Coal feed rate (number 
basis) multiplied by 
and divided by char withdrawal 
rate (number basis) 

^cw 

BC 


Dimensionless ratio (see 
ref. 4) 


BETA 


Adjustable parameter for 
elutriation calculations 



BEDVOL 

Total bed volume 

c 



Wet coal feed rate divided by 
air feed rate 


CABS 

CABS 

Mole ratio of calcium to 
sulfur in the feed solids 



CADF 

Heat capacity of limestone fed 



CCF 

Heat capacity of coal fed 


CELU 


Elutriation rate of solid 
carbon 
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Mathematical 

Level I 
FORTRAN 

Level II 
FORTRAN 




CGMF 

Molar heat capacity of feed 
gas at feed temperature 



CLOSS 

Total carbon loss to solids 

Db 

DB 


Average bubble diameter 

DbM 



Maximum bubble diameter 

DbO 



Initial bubble diameter 


DPL 


Mean diameter of limestone 
particles 




Mean diameter of limestone 
particles in size interval i 


DPF 

DPCF 

Diameter of coal fed to 
combustor at various sieve 
sizes 


DT 


Bed diameter 



DNZL 

Diameter of distributor holes 



DPAOF 

Diameter of limestone fed to 
combustor at various sieve 
sizes 



DTHICK 

Distributor plate thickness 



DTUBE 

Outside diameter of heat 
exchanger tubes 



DZAV 

Specified compartment height 
for calculations 

E 

EXAIR 

EXAIR 

Excess air ratio 


EC 


Quantity used in elutriation 
calculations (see ref. 4) 


ELOSS 

ELLOSS 

Fraction of burnable carbon 
elutriated 


EMF 


Void fraction of bed at 
minimum fluidization 
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Mathematical 

Level I 
FORTRAN 

Level II 
FORTRAN 




ETS 

Efficiency of sulfur capture 

^m 

FMO 


Total molar flow rate of gas 
in the bed 



FW 

Volume ratio of wake to bubble 


FRACT(I) 

FRACTC(I) 

Weight fraction of coal fed to 
combustor between DPF(I-l) and 
DPF(I) 



FD 

Fraction of solids discharged 
from bed at each discharge 
location ^excluding elutriated 
solids) 



FFAD 

Fraction of limestone fed at 
each solids feed location 



FFC 

Fraction of coal fed at each 
solids feed location 



FRACTA(I) 

Weight fraction of limestone 
fed to combustor between 
DPADF(I-l) and DPADF(I). 



FS 

Fractional conversion of 
1 imestone 

g 



Acceleration of gravity 



GZCO 

Ratio of moles of CO formed to 
moles of C in coal fed 



6ZH2 

Ratio of moles of H 2 formed 
to sum of moles of H 2 
formed, moles of H 2 O formed, 
and moles of moisture in coal 
fed 



GZH2S 

Ratio of moles of H 2 S formed 
to total moles of S in feed 


HCHAR 


Hold-up of char in bed 


HRC 


Mass fraction of carbon in bed 
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Mathematical 

Level I 
FORTRAN 

Level II 
FORTRAN 




HAREA 

Total heat transfer area 



HLF 

Height of the fluidized 
(expanded) bed 


lELUTR 


Selector for elutriation 
correlation (1 for ref. 15, 
2 for ref. 16, or 3 for 
ref. 17) 


IMODEL 


Selector for gas flow model 
(2 for plug flow) 



lARR 

Selector for arrangement of 
heat exchanger tubes (0 for 
none, 3 for horizontal inline) 



IGNITE 

Selector for cold flow or 
combustor (1 for combustion) 

*-mf 

HLMF 

HLMF 

Minimum fluidization bed 
height 

Mair 



Molecular weight of air 

Me. Mca. 
Mh2> Mh20. 
Mn2> Mq2. 
Ms 



Molecular weight of species 
indicated 



MDIS 

Number of solids discharge 
locations 



MFEED 

Number of solids feed 
locations 



MTB 

Number of heights where 
ATB is given 



MTHE 

Number of heat exchanger 
sections, including sections 
with none 

nd 


AND 

Number of distributor holes 


N 

NDPC 

Number of members of FRACT 
or FRACTC array 



Mathematical 


Level I 
FORTRAN 


Level II 
FORTRAN 
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NDPAD 

Number of members of FRACTA 
array 

p 

P 

PAV 

Average absolute pressure in 
bed 


PASH 


Fraction of ash elutriated 



PF 

Absolute gas pressure below 
the distributor 



PH 

Horizontal pitch of heat 
exchanger tubes 



PV 

Vertical pitch of heat 
exchanger tubes 



QAREA 

Heat transfer rate to 
heat exchangers per unit 
area 



QCOAL 

Lower heating value of wet 
coal 



QTRANS 

Total heat transferred to 
heat exchangers 



QVOL 

Heat transfer rate to heat 
exchangers per unit volume 
of bed 

R 



Universal gas constant 


RHOC 


True density of coal 



SOLVOL 

Apparent volume of solids 
in bed 

T 

TK 

TAV 

Average bed temperature 
(initial estimate in case 
of Level II) 



T 

Bed temperature 



TETUBE 

Total volume fraction of heat 


exchanger tubes 
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Mathematical 

Level I 
FORTRAN 

Level II 
FORTRAN 




TF 

Inlet gas temperature 



TWAV 

Average heat exchanger coolant 
temperature 




Superficial gas velocity at 
minimum fluidization 

Uo 

UO 


Superficial gas velocity in 
bed at T and p 



UHEAV 

Average overall heat transfer 
coefficient between bed and 
heat exchanger coolant 



UF 

Superficial gas velocity at 
the distributor at TF and PF 

'^cone 



Volume of fluidized conical 
bed 


VCF 


Volumetric flow of coal feed 
excluding voids 


VCW 


Volumetric withdrawal rate of 
solids excluding voids 


VLS 


Volumetric feed rate of 
limestone excluding voids 


VMF 

VMF 

Volume of bed at minimum 
fluidization everywhere' 
(hypothetical for conical 
bed) 

Wc 

WCOAL 


Dry coal feed rate by weight 


WBED 


Total bed weight 



WAD 

Total limestone feed rate by 
weight 



WCOAL 

Wet coal feed rate by weight 



WDIS 

Total discharge rate by weight 



Mathematical 


Level I 
FORTRAN 


Level II 
FORTRAN 
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WELT 

Solids elutriation rate by 
wei ght 

Xi 



Weight fraction in size 
interval i 

K 



Weight fraction total carbon 
in dry coal 

XH 

XH 

XH 

Weight fraction hydrogen in 
dry coal 

V,Ca 

XLCA 


Weight fraction calcium in 
1 imestone 

Xm 

XW 

XW 

Weight fraction moisture in 
coal (dry basis) 

XN 

XN 

XN 

Weight fraction nitrogen in 
dry coal 

XO 

XO 

XO 

Weight fraction oxygen in dry 
coal 

XS 

XS 

XS 

Weight fraction sulfur in dry 
coal 


XCF 

XCF 

Weight fraction fixed carbon 
in dry coal 


XCV 

XCV 

Weight fraction volatile 
carbon in dry coal 


XGO 


Mole fraction oxygen in 
effluent gas 



XAV 

Average burnable carbon 
concentration in the bed, 
weight basis 



XCAC03 

Weight fraction CaC 03 in 
1 imestone 



XGF 

Feed gas composition 



XGO(l) 

Mole fraction O 2 in 


effluent gas 
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Mathematical 

Level I 
FORTRAN 

Level II 
FORTRAN 




XGO(2) 

Mole fraction CO 2 in 
effluent gas 



XGO{3) 

Mole fraction SO 2 in 
effluent gas 



XGO(4) 

Mole fraction H 2 O in 
effluent gas 



XM6C03 

Weight fraction MgC 03 fn 
1 imestone 

y 



Ratio of char diameter to 
maximum char diameter 

yo2,o 

XGOO 


Mole fraction O 2 in the feed 
gas 

^cyl 



Bed depth of equivalent 
cylindrical bed 



ZB 

Height above distributor 
where ATB specified 



ZDIS 

Height of solids discharge 
above distributor 



ZF 

Height of solids feed above 
distributor 



ZHE 

Height of tops of heat 
exchanger sections (including 
sections with no heat 
exchangers) 




Volume fraction of bubbles 


ETC 

ETC 

Combustion efficiency 



ETCA 

Initial estimate of 
combustion efficiency 

(9 

THET 


Mean residence time of 
limestone 
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Mathematical 

Level I 
FORTRAN 

Level II 
FORTRAN 





Viscosity of gas in bed 

pg 



Density of gas in bed 

P£S 

RHOLS 

RHOAO 

True density of limestone 




Bed particle density 




Distribution function for 
char in bed, number basis 


Note: "Wet" coal means coal as received from producer. 
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APPENDIX B - METHODS OF CALCULATING OR OBTAINING 
INPUTS TO PROGRAMS FOR CASES la, lb, Ic, 2a, 3a, AND 4a 
This appendix gives methods of calculating or obtaining inputs to 
Level I and II for the five averaged tests for those inputs which were 
the same for Cases la, lb, Ic, 2a, 3a, and 4a. Values are also given. 

Level I Inputs 
XCF 

This was obtained from the proximate analysis of coal and was 
0.5340. 

XCV 

This was obtained from the proximate and ultimate analyses of coal 
and was 0.2198. 

XH, XS. XO, and XN 

These were obtained from the ultimate analysis of coal and were 
0.0514, 0.0199, 0.0761, and 0.0149, respectively. 

xw 

This was obtained from the proximate analysis of coal and was 
0.0217. 

N, DPF. and FRACT 

These were obtained by dry sieving coal taken from between the 
metering screw and the blending auger. N was 19, and DPF as well as 
FRACT are given in Table VII. 
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WCOAL 

This was obtained from the rotation rate of the fuel metering screw 
and the limestone-to-coal ratio and was 4.51 g/sec. 

CABS 

The limestone-to-coal ratio was obtained from rotation rates of the 
limestone and coal metering screws. Then CABS was calculated from the 
limestone-to-coal ratio and chemical analyses of the limestone and coal 
and was 1.617. 

EXAIR 


This was obtained from the equation (ref. 4) 



and had the value 0.639. 
VO 


This was 0 since EXAIR was specified. 

P and TK 

These were observed directly and had values of 5.15 atm and 1151K, 
respectively. 

This was taken from a drawing and was 36.4 cm. 

PASH 

This was measured with help from chemical analysis of fly ash and 
discharge solids and was 0.780. 
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IMODEL 

This was 2 for plug flow of gas. 

EMF 

This was measured on spent bed material from the combustor in the 
cylindrical bed (fig. 2) and was 0.51. 

RHOC 

This was measured on the coal in a pyknometer with Jet A kerosene 

3 

and was 1.267 g/cm . 

XLCA 

This was obtained from chemical analysis of the limestone and was 
0.383. 

Level II Inputs 

In addition to the quantities already given for Level I, some of 
which are also used in Level II: 

CADE 

From reference 20 the specific heat of limestone is 0.217 cal/g C. 
CCF 

From references 21 and 22 using the Kirov approximation, the 
specific heat of the wet (as received) coal is 0.3139 cal/g K. 

CGMF 

From references 23 and 24, the specific heat of air is 6.955 
cal/mole K. 

MTB, ZB, and ATB 

From drawings MTB is 4, ZB(1) - ZB(4) are 0., 62.7, 81.3, and 280. 
cm; and ATB(l) - ATB(4) are 405., 405., 670., and 2193 cm^. 
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MTHE, ZHE, AHE. DTUBE. PV, PH, and lARR 

MTHE was 5. See Table VIII for values of other variables from 
drawings. 

MFEED. ZF(1). FFC(l), and FFAD(l) 

From drawings the values were 1, 6.91 cm, 1., and 1., respectively. 
MDIS. ZDIS(l), and FD(1) 

From drawings the values were 1, 141.94 cm, and 1., respectively. 
AND. DNZL, and DTHICK 

The values were 36., 0.317 cm, and 1.27 cm, respectively. 

DZAV 

To be consistent with Table VIII, a value of 8 cm was used. 

FW 

The value of 0.15 in the test case (ref. 4) was used. 

XCAC03 and XMGC03 

Values were obtained from chemical analysis of the limestone and 
were 0.9575 and 0.0107. 

QCOAL 

The higher heating value was measured in an ASTM bomb calorimeter 
and converted to a lower heating value of 7109 cal/g. 

NDPC. DPCF. and FRACTC 

NDPC was 19. DPCF and FRACTC are given in Table VII. 

HLF, VMF. and HLMF 

HLF from drawings was 141.94 cm (only one of these three quantities 
need be given; the others are then set to 0). 
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PAV and TAV 

Same as P and TK in Level I. 

TWAV 

Measured as 298K. 

WCOAL 

Calculated the same as in Level I except for wet (as received) 
basis, so value was 4.62 g/sec. 

WAD and CABS 

CABS was same as in Level I (only one of these two quantities need 
be given; the other is then 0). 

This is always set to 0 if EXAIR is specified. 

TF and PF 

Measured as 297K and 5.47 atm, respectively. 

XGF 

For air (ref. 24) the values of XGF(1)-XGF(7) were 0.21, 0., 0., 
0., 0., 0., and 0., respectively. 

6ZC0 and GZH2 

These were measured as 4.45x10"^ and 0., respectively. 

6ZH2S 

Since all cases were lean, this was assumed to be 0. 

IGNITE 

Assuming combustion, the value was 1. 
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APPENDIX C - METHODS OF CALCULATING OR 
OBTAINING INPUTS TO PROGRAMS 
FOR CASES la, lb. AND Ic 

This appendix gives methods of calculating or obtaining inputs to 
Level I and II for the five averaged tests for those inputs which were 
not the same for Cases la, lb, Ic, 2a, 3a, and 4a and applies to Cases 
la, lb, and Ic only. Values are also given. 

Level I Inputs 
PPL 

The limestone was dry sieved (see Table VII). The mean diameter 
(DPL) was then found from (ref. 4) 

(Cl) 

and was 0.0790 cm. 

RHOLS 

The true density of the limestone was measured in a pyknometer with 

•5 

Jet A kerosene and was 2.66 g/cm . 

First F^ was found from the equation (ref. 4). 
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Next was found from 




F 

fpl 


(C3) 


Viscosity was found from (ref. 4) 


/-C = 3.72 X 10"® T 


fcr) 


Gas density was found from 


f: 


pM . 


/ FT 


(C5) 


Bed particle density was found from (ref. 4 with change from Table I) 




X 




(C6) 


^JS 




LA 


hr: 


/ = 


Minimum fluidization velocity was found from (ref. 4) 

— ^ I 33,7 t . 






7>(C7) 


Maximum bubble diameter was found from (ref. 4) 


D 


Bm 


= 0.652[a^(u„-u„,^'>-' 


Initial bubble diameter was found from (ref. 4) 

D zo377\—^ — ^ 

n . 


(CF) 


(C9) 


ORIGINAL 

m 
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The bed depth for the equivalent cylindrical bed was found from 





rd/1 e. 


(ClOl 


Finally Dg was found from (ref. 4) 

~D 

D - D, - ( D - A- £ ^ 

and was 16.00 cm. 

HLMF 

First the bubble fraction was found from (ref. 4) 



I 

I ^ 

~ u r 
o f^T 


(Cll) 


(C12) 


Then the minimum fluidization bed height (HLMF) was found from (ref. 4) 


L 


V 

C Ol^t 




(C13) 


A- i 


and was 57.5 cm. 

BETA 

This is a parameter to adjust to match the experiment. It was 
taken to be one because this would be its value if the elutriation 
correlations were directly applicable. 
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Level II Inputs 

RHOAD 

Same as RHOLS in Level I inputs. 

NDPAD. DPADF. and FRACTA 

NDPAD was 19. DPADF and FRACTA are given in columns 1 and 3 of 
Table VII. 

UHEAV 

-3 7 

Measured as 7.65x10 cal/sec cm K. 
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TABLE I. - CHANGES TO LEVEL I PROGRAM 


Program or 
subprogram 

Purpose 

Change 

ALL 

Obtain sufficient numerical 
accuracy so 9 will converge 

Run on UNIVAC 1100 series computer 
instead of IBM 360 (both in single 
precision) . 

MAIN 

Obtain consistent value of 
density of calcined limestone 

Change card 22 from RHOL = RHOLCA to 
RHOL = RHOLCA * 56.08 / MCA 

MAIN 

Eliminate use of 0 for two 
different quantites at same i 

time 

Change card 57 from 
D = 4.26 * (TK / 1800.) ** 1.75 / P to 
002 = 4.26 * (TK / 1800.) ** 1.75 / P 

MAIN 

Put particle diameters in 
increasing order to facilitate | 
numerical integration 1 

Delete card 128. 

Change card 130 from IF(K.EQ.N) GO TO 21 
to 16 IF(K.EQ.N) GO TO 21 
Insert after card 132 

IF(XX.GT.DPF(K)) GO TO 19 
DP(L) = XX 

Insert after card 133 
GO TO 17 

19 DP(L) = DPF(K) 

PHIF(L) = FRACT(K) 

GO TO 16 

MAIN 

Eliminate use of D for two 
different quantities at same 
time 

Insert after card 142 
D = 002 

MAIN 

Obtain correct mass feed 
rate of uncalcined lime- 
stone 

Change card 201 from 

WLS = CABS * WCOAL * XS * MCA / MS to 

WLS = CABS * WCOAL * XS * MCA / (MS * XLCA) 

MAIN 

Obtain correct volumetric 
feed rate of limestone 

Change card 202 from 
VLS = WLS / RHOLCA to 
VLS = WLS / RHOLS 

MAIN 

Obtain e convergence 

Change card 211 from 
DTH = -THETM * 0.25 to 
DTH = -THETM * 0.1 

MAIN 

Obtain 9 convergence 

j 

Change card 225 from 

THET = THETM / (1 . + THl * ETC) to 

THET = THETM 

MAIN 

Obtain convergence 

Change card 243 from 
DBCl = 100. to 
DBCl = 10. 

Change card 244 from 
DO 50 IBCl = 1, 20 to 
DO 50 IBCl = 1, 200 

MAIN 

Eliminate use of A1 for two 
different quantities at same 
time 

Change card 288 from 
A1 = AT / (WBED * EC) to 
AAAA = AT / (WBED * EC) 

Change card 301 from 

YCl = DP(I) ** 3 * A1 * AKE(I) * PHI(I) to 
YCl = DP(I) ** 3 * AAAA * AKE(I) * PHI(I) 

POP 

Eliminate negative values 

of < 1 ) 

Insert after card 29 

IF(PHI(I).LT.O. ) PHI(I) = 0. 


original L \GE 

OE POOR QUALIXa 







TABLE II. - CHANGES TO LEVEL II PROGRAM 


Program or 
subprogram 

Purpose 

Change 

Main 

Make N0« correlation agree 

Change card 468 from 


with reference 14 

ANOX = AN * WCOAL * (1. - ETC) * QCOAL / FMO to 
ANOX = AN * WCOAL * ETC * QCOAL / FMO 

HYDRO 

Obtain correct values of bed 

Insert after card 843 


volume 

BEDVOL = BEDVOL - DVBB(I) 

VOLUME 

Get HLMF, VMF, and HLF to be 

Change card 954 from 


consistent no matter which is 

A = 1.0 - (ZZ - Z(N)) / DZAV to 


nonzero on input 

A = - (ZZ - Z(N)) / DZAV 







TABLE III. - OUTPUT FOR LEVEL I TEST CASE 


Quantity 

(see 

appendix A) 

Value from 
reference 4 
or IBM 360 
computer at NASA 

Value with 
changes in table I 
using UNIVAC 1100 
series computer 

ETC 

0.9291 

0.9327 

XGO 

0.5332x10"^ 

0.3316 

THET, sec 

0.1859x10^ 

0.1793x10® 

BC 

0.1005 

,0.1070 

BCl 

0.1701x10^ 

0.1650x10® 

VCW, cm^/sec 

0.8210x10'^ 

0.1889 

VCF, cm^/sec 

0.2086x10^ 

0.2086x10® 

VLS, cm^/sec 

0.2309x10^ 

0.2309x10^ 

HCHAR, g 

0.1386x10^ 

0.3077x10^ 

UBED, g 

0.7832x10® 

0.9296x10® 

HRC 

0.1375x10'^ 

0.2571x10'^ 

EC, 1/sec 

0.1131x10'® 

0.1065x10'® 

ELOSS 

0.6694x10'^ 

0.5821x10'^ 

ALAMV 

0.1236x10'® 

0.2971x10'® 

ALAMA 

0.7792x10'® 

0.1081x10'® 








TABLE IV. - OUTPUT FOR LEVEL II TEST CASE 


— 

Quantity 
(see appendix A) 

Value from 
reference 4 

Value with 
changes in table II 
using UNI VAC 1100 
series computer 

ETC 

0.9282 

0.9282 

XAV 

0.1547x10’^ 

0.1547x10'^ 

TAV, K 

0.7759x10^ 

0.7757x10^ 

ETS 

0.8287 

0.8288 

FS 

0.3767 

0.3767 

XGO(l) 

0.4560x10'^ 

0.4560x10'^ 

XG0{2) 

0.1635 

0.1635 

XGO(3) 

0.3500x10’^ 

0.3498x10'® 

XG0(4) 

0.5742x10'^ 

0.5742x10'^ 

ANOX 

0.1839x10"'^ 

0.2377x10'® 

HLF, cm 

0.6705600x10^ 

0.6706x10® 

HLMF, cm 

0.4842171x10^ 

0.4843x10® 

VMF, cm^ 

0.1791724x10® 

0.1792x10® 

BEDVOL. cm^ 

0.3013810x10® 

0.3014x10® 

SOLVOL, cm^ 

0.1791724x10® 

0.1792x10® 

TETUBE 

0.1131662 

0.1132 

HAREA, cm^ 

0.2526380x10® 

0.2526x10® 

QTRANS, cal/sec 

0.7488019x10® 

0.7485x10® 

QVOL, cal/sec cm^ 

0.2484568 

0.2484 

QAREA, cal/sec cm^ 

0.2963932x10^ 

0.2963x10^ 

WELT, g/sec 

0.2945155x10^ 

0.2945x10^ 

CELU, g/sec 

0.7362888 

0.7363 

CLOSS, g/sec 

0.7550509 

0.7551 

WDIS, g/sec 

0.1215159x10^ 

0.1215x10^ 








TABLE V. - IMPORTANT INPUT AND OUTPUT FOR LEVEL I PROGRAM FOR PRELIMINARY RESULTS (CASE la, lb & Ic) AND CHANGES (CASES 2a, 3a & 4a). 


OBSERVED VALUES FOR THE FIVE AVERAGED TESTS ARE GIVEN FOR COMPARISON 


Case 

Inputs 

Outputs 


Minimum 

Mean 

Adjustable 

Average 

Limestone 

Selector 

Combustion 

Limestone 

B X coal 
cw 

Fraction of 


fluidization 

diameter 

parameter 

bubble 

true 

for 

efficiency. 

residence 

feed rate 

burnable 


bed 

of 

for 

diameter. 

density. 

elutriation 

ETC 

time. 

T char 

carbon 


height, 

1 imestone. 

elutriation. 

DB 

RHOLS 

correlation, 


THET 

withdrawal 

elutriated, 


HLMF 

DPL 

BETA 

cm 

g/cm^ 

lELUTR 


sec 

rate. 

ELOSS 


cm 

cm 







BCl 


la 

57.5 

0.0790 

1 . 


2.66 


0.999108 

1.231x10^ 

■n 

8.08x10"^ 

lb 

57.5 

,0790 

1. 


2.66 

‘’2 

.999846 

1.231x10^ 


6.15x10’® 

Ic 

57.5 

.0790 

1. 

■■■ 

2.66 


.999885 

1 .231x10^ 

msm 

1.57x10’® 

2a 

62.4 

.0793 

1 . 


4.89 

"l 

.999455 

2.03x10^ 

13.10 

4.91x10’^ 

3a&4a 

62.4 

.0793 

25. 

Hm 

4.89 

"l 

.990778 

2.04x10^ 

17.36 

1.048x10’^ 

Observed 


.0790 


■H 

2.66 


.9895 



1.048x10’^ 


^Elutriation per reference 15. 
^Elutriation per reference 16. 

r 

Elutriation per reference 17. 























TABLE VI. - IMPORTANT INPUT AND OUTPUT FOR LEVEL II PROGRAM FOR PRELIMINARY RESULTS (CASE la) AND CHANGES (CASES 2a, 3a & 4a). 


OBSERVED VALUES FOR THE FIVE AVERAGED TESTS ARE GIVEN FOR COMPARISON. 


Case 


Inputs 




Outputs 




Limestone 

Size 

Overall 

Combustion 

Bed 

Efficiency 

Fractional 

Mole 

Mole 


true 

distribution 

heat 

efficiency. 

tempera- 

of sulfur 

conversion 

fraction 

fraction 


density, 

RHOAD 

g/cm^ 

of 

limestone, 
NDPAO, 
DPADF, 
FRACTA for 

transfer 

coefficient, 

UHEAV 

2 

cal /cm 
sec K 

ETC 

ture, 

t" 

C 

capture, 

ETS 

of 

limestone, 

FS 

SO 2 in 
effluent, 
XG0(3) 

NO^ in 
effluent, 
ANOX 

la 

2.66 

limestone 

7.65x10'^ 

0.9990 

980 

0.611 

0.378 

0.000408 

0.000290 

2a 

4.57 

spent bed 

7.65x10'^ 

0.9994 

1005 

.785 

.485 

.000225 

.000290 

3a 

4.57 

spent bed 

7.65x10'^ 

0.9876 

1019 

.834 

.515 

.000175 

.000287 

4a 

4.57 

spent bed 

1.083x10"^ 

.9876 

878 

.861 

.532 

.000146 

.000287 

Observed 

2.66 

1 imestone 

7.65x10'^ 

.9895 

878 

.792 

.490 

.000226 

.000125 


^Temperature at 36 cm above top of distributor. 
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TABLE VII. - RESULTS OF DRY SIEVING COAL. LIMESTONE, AND SPENT BED 


MATERIAL. QUANTITIES ARE INPUTS TO LEVEL I AND II PROGRAMS 


Diameter 
at top of 
size 

interval , 
DPF, 
DPADF, 
or DPCF for 
all cases, cm 

Weight 

fraction coal 
in size 
interval , 
FRACT or 
FRACTC for 
all cases 

Weight 

fraction lime- 
stone in 
size interval, 
FRACTA® 
for cases la, 
lb & Ic 

Weight 

fraction spent 
bed in size 
interval , 
FRACTA*^ 
for cases 2a, 
3a, 4a & 
5-30 

0.0074 

0.0003 

0.0123 


.0104 

.0009 

.0068 


.0147 

.0069 

.0135 


.0175 

.0025 

.0026 


.0208 

.0163 

.0032 


.0295 

.0356 

.0088 


.0351 

.0425 

.0053 

0.0170 

.0417 

.0538 

.0091 

.0239 

.0495 

.0553 

.0114 

.0508 

.0589 

.0594 

.0123 

.0645 

.0701 

.0825 

.0255 

.1515 

.0833 

.0766 

.0622 

.1645 

.0991 

.1172 

.1424 

.1998 

.1168 

.0963 

.1259 

.1246 

.1397 

.0960 

.1406 

.0938 

.1651 

.1500 

.1856 

.0827 

.1981 

.0747 

.1433 

.0213 

.2362 

.0194 

.0854 

.0049 

.2794 

.0138 

.0038 

.0007 


®For limestone from between limestone metering screw and 
blending auger. 

*^For spend bed material. 




TABLE VIII. - SPECIFICATIONS OF HEAT EXCHANGER TUBES FOR LEVEL II 


PROGRAM INPUT FOR ALL CASES 
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Figure 1. - Schematic of leRC pressurized fluidized bed combustor. 











Figure 3. - Comparison of burnable carbon elutriated based on experiment and 
three eiutriation correiations for preliminary resuits (BETA 1). 









Figure 8. - Comparison of observed and calculated nitrogen oxides emission rate 
tor cases 5 to 30. 



Figure 9. - Comparison of observed and calculated bed temperature at 14. 7 
inches above the distributor for cases 5 to 30. 



0 .2 .4 .6 .8 1.0 1.2 1.4 

Observed sulfur emission rate, lb 502/10^ Btu 

Figure 10. - Comparison of observed and calculated sulfur emission rate for 


cases 5 to 30. 



Observed efficiency of sulfur capture 


Figure 11. - Comparison of observed and calculated effi- 
ciency of suifur capture for cases 5 to 30. 
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Figure 13. - Effect of bed depth on ratio of caiculated to observed 
burnable carbon eiutriated for cases 5 to 30 (theoretical elu- 
triation (ref. 15) multiplied by BETA - 25). 


